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We analyze Tevatron data of bottomonium hadroproduction in the framework of the colour-octet model (COM) 
implemented in the event generator PYTHIA |lj using CTEQ4L PDF taking into account initial-state radiation 
of gluons and Altarelli-Parisi evolution of final-state gluons. We obtain new values for the colour-octet matrix 
elements (Me's) relevant to this production process for the T(nS) family (n=l,2,3), finding that 1 Sq + 3 Pj 8 ' 
contributions are not needed in the fit. We show the different contributions to Y(IS') production at Tevatron for 
Pt > 8 GeV, comparing them with CDF data. Remarkably we find a quite small contribution (compatible with 
zero) from feeddown of Xbj states produced through the colour-octet mechanism: T(IS') indirect production via 
Xbj decays should be mainly ascribed to the colour-singlet model. Finally we extrapolate to LHC energies to 
predict T(nS) production rates. 



1. HIGHER ORDER QCD EFFECTS; 
GENERATING T(nS) 

In a previous work we extracted colour- 
octet Me's for T(IS') hadroproduction using 
CTEQ2L parton distribution funcion (PDF) @; 
in this work we extend this study to the upper 
resonances T(2S) and T(3S*), using the improved 
CTEQ4L PDF {§. We have implemented COM 
in the event generator PYTHIA with initial-state 
radiation of gluons |5|,|6| and Altarelli-Parisi evo- 
lution of final-state gluons in the same way that 
earlier. 

First we repeat our fit of the Y(IS) differen- 
tial cross section with CTEQ4L PDF. We gener- 
ated separately each source: direct T(15) produc- 
tion from both Colour Singlet Model (CSM) and 
COM. Indirect source comes from electromag- 
netic decays of x&j(1P) and Xbj(2P) and contri- 
butions from strong and electromagnetic decays 

of T(2S) and T(3S); the M 5 = 5x ( < ° s(3 f ,)> + 
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<o 8 (^ g o)> ^ combination is also taken into ac- 
count (although later we shall see that actually 
this contribution is not needed). For the CSM 
Me's we take the values from the Buchmuller and 
Tye QCD potencial jjj. We make use of the 
Heavy Quark Spin Symmetry in order to relate 
the different matrix elements. 

We set the T(15) mass as 9.46 GeV (the mass 
of the coloured intermediate state was set as 
2Mb = 9.76 GeV, the difference is accounted for 
by soft gluon emission) whereas for Xbj{ n P) we 
took a weighted mean of their real masses || . For 
T(25) we consider direct CSM contribution and 
the decay from T(35) and Xbj(2P)- COM con- 
tribution also is taken into account. For Xb,j(2P) 
we set its real mass. Regarding the T(35) state 
we consider only direct production for both CSM 
and COM contributions; we set its real mass, ex- 
cluding any Xbj(3P) decays. 

We have included the higher order QCD effects 
on the partonic cross sections caused by: 

• Intrinsic Fermi motion of partons inside the 
hadrons. This non perturbative effect is rel- 
evant at small P t values. 



• Multiple emission of gluons in the initial 
state. This perturbative contribution is dy- 
namically generated via gluon radiation im- 
plemented in the event generator PYTHIA 
by means of a parton shower algorithm. 
This effect overshadows the former at high 
Pt values. 

• Altarelli-Parisi (AP) evolution of the split- 
ting gluon in the gg^bbg channel. Actually 
we do not generate the virtual gluon (<?*) 
that splits into bb pairs; final hadronization 
into a (QQ) bound state is taken into ac- 
count by means of the colour-octet matrix 
elements multiplying the respective short- 
distance cross sections. Nevertheless, it is 
reasonable to assume that, on the average, 
the virtual g* should evolve at high pt sim- 
ilarly to the other final-state gluon - which 
actually is evolved by the PYTHIA machin- 
ery. We used this fact to simulate the (ex- 
pected) evolution of the (ungenerated) g* 
whose momentum was assumed to coincide 
with that of the resonance (neglecting the 
effect of emission/absorption of soft gluons 
by the intermediate coloured state bleeding 
off colour [§). 

2. RESULTS: 

In table 1 we present the results of our fit to 
CDF data ||. Table 2 shows our previous result 
using CTEQ2L (both in 1CT 3 GeV 3 units); errors 
are statistical only. Figures 1 and 2 show the the- 
oretical curves obtained for T(1S) and T(3S) re- 
spectively. The statistical x 2 /Ndf are quite good, 
especially for T(35). Those values are consistent 
with NRQCD velocity scaling rules. Let us re- 
mark that due to the pt cut-off parameter set in 
the generation, only those experimental points for 
Pt > 1 GeV were used in the fit. 

Tablet. COM parameters (CTEQ4L): 



<0 8 T(1S) ( 3 S!)> 


\tot 


= 80 ± 28 


<0 8 T(2S) ( 3 5 1 )> 


\tot 


= 77 ± 42 


<0 8 T(3S) ( 3 Si)> 


tot 


= 90 ± 32 



Table2. COM parameters (CTEQ2L): 

<0 8 T(1S) ( 3 5 1 )>| fot = 139 ±31 
<0 8 T(2S) ( 3 5 1 )>| tot = 80 ±44 
<0 8 T(3S) ( 3 5 1 )>| tot = 75 ±27 
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Figure 1. Theoretical curve obtained from the fit 
using CTEQ4L PDF in T(15) with AP evolution. 



As a check of our generation we reproduce the 
relative fractions of the different contributions to 
T(15) production obtained by the CDF Collabo- 
ration |0. In table 3 we quote the results from 
CDF data at P T > 8 GeV; in table 4 the re- 
sults from our fit in the same case. The agree- 
ment is quite good taking into account the errors. 
If we compare experimental data versus genera- 
tion, we can see that the 3 s[ 8S> contribution goes 
mainly to direct T(IS') production and therefore 
the Xbj( n P) channel would be mainly saturated 
by the CSM. As a consequence the T(15) COM 
matrix element should be mainly direct. 
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Figure 2. Theoretical curve obtained from the fit 
using CTEQ4L PDF in T(35) with AP evolution. 



Fitting the Y(35) differential distribution we 
assumed only direct production (i.e. no feeddown 
from higher resonances) and hence the COM pa- 
rameter must be entirely ascribed to direct con- 
tribution. If we compare with the previous val- 
ues obtained by Cho and Leibovich for direct 
production: 5.9; 4.1; 4.1 (in 10~ 3 GeV 3 units) for 
T(nS) (n = 1,2,3 respectively) we find them too 
low, whereas their corresponding values (consid- 
ering now all contributions) < Os( 3 Si) > \ tot too 
high: 480; 220; 160 ( in 10~ 3 GeV 3 units) respec- 
tively, since the COM Xbj{nP) contributions are 
overestimated according to our analysis. 

Table 3. Rel. fract.(in %) for T(15) from CDF 

T(lS)L rect = 50.9 ±12.2 

T(25) + T(35) = 11.5 ±9.1 

XbjilP) = 27.1 ±8.2 

Xb.j(2P) = 10.5 ±4.6 



Table 4. Rel.fract.(in %) for T (15) (generation) 

T(15)| 3S (8, +CSM = 32.7 + 20.7 -52.4 

T(25)+T(35)| CSM = 4.1 

Xbj(lP)\ CSM = 25.7 

Xbj(2P)\ CSM = 16-8 



3. EXTRAPOLATION TO LHC 

We predict the differential and integrated cross 
section at the LHC taking the < Os( 3 5i) > ma- 
trix elements obtained from our fit to CDF data. 
We have found gg — gq — qq contributions at LHC 
energies to be 77%, 22%, 1% respectively. Figures 
3 and 4 show the differential cross section at LHC 
energies for T(15) and T(35) respectively. We 
can see that T(35) should be the best candidate 
to check the COM since, even at relatively low 
Pt, the 3 s[ 8 ^ contribution becomes dominant. In 
table 5 we summarize the results of the integrated 
cross section for each Y(n5) (IS total and direct 
production). 

Table 5. Integrated cross sections (nb) 



T(15)| tot 


<T = 


33.4 


ms)\ dir 


a — 


12.4 


Y(25)| tot 


a = 


9.4 


Y(35)| iot 


a = 


4.6 



4. CONCLUSIONS 

• We have analyzed T(n5) hadroproduction 
at the Tevatron concluding that: 
< 0^ nS \ 3 Si) > slightly varies in chang- 
ing from CTEQ2L to CTEQ4L ( as ex- 
pected since both PDF's have a different 
behaviour, especially at low x) and are 
compatible with NRQCD velocity scaling 
rules. We find the Xbj(nP)\ COM contribu- 
tion very small. Thereby comparing with 
Cho and Leibovich work, all < 0$( 3 Si) > 
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Figure 3. Theoretical curve using CTEQ4L PDF 
in Y(IS') with AP evolution. 



Figure 4. Theoretical curve using CTEQ4L PDF 
in T(3S) with AP evolution. 



values for direct production have to be en- 
hanced. 

. M 5 = 5x ( <0s ff j)> + <O8( ; S0)> ) is com- 
patible with zero. 

• We have extrapolated to LHC energy pro- 
viding theoretical curves for the differential 
transverse momentum distribution and the 
total cross section for each Y(nS). Produc- 
tion rates from low to high P t will afford 
very interesting physics for bottomonium 
at the LHC, especially polarisation for the 
T(35) resonance at p t > 20 GeV. 
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